Measurement of interaction energy near a Feshbach resonance in a 6 Li Fermi gas 
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We investigate the strongly interacting regime in an optically trapped 6 Li Fermi mixture near 
a Feshbach resonance. The resonance is found at 800(40) G in good agreement with theory. 
Anisotropic expansion of the gas is interpreted by collisional hydrodynamics. We observe an unex- 
pected and large shift (80 G) between the resonance peak and both the maximum of atom loss and 
the change of sign of the interaction energy. 
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The achievement of Bose-Einstein condensation (BEC) 
in dilute atomic gases Q has naturally triggered research 
on cooling of Fermi gases to quantum degeneracy. Sev- 
eral groups have now reached Fermi degeneracy in 40 K 
and 6 Li with temperatures down to about 0.1 to 0.2 of 
the Fermi temperature T F 0, 0, 0, 0, H, . One of the 
major goals of this research is to observe the transition 
to a superfiuid phase the analog of the supercon- 
ducting (BCS) phase transition in metals |9j. Very low 
temperatures and strong attractive interactions in a two- 
component Fermi gas are favorable conditions to reach 
this superfiuid state. The interactions in atomic gases at 
low temperature are usually described by a single param- 
eter, the s-wave scattering length a. This quantity can 
be tuned near a Feshbach resonance where the sign and 
magnitude of a can be adjusted by means of an exter- 
nal magnetic field B |10L lll|. Therefore, when entering 
the regime of strong interactions, dilute Fermi gases offer 
unique opportunities to test new theoretical approaches. 
For instance, near resonance, the critical temperature for 
superfluidity has been predicted to be as high as 0.25-0.5 
Tp |l2^. a temperature range experimentally accessible. 

Both 40 K and 6 Li possess Feshbach resonances at 
convenient magnetic fields [H El In Fig.Cl is plot- 
ted the theoretical scattering length for the mixture of 
the two lowest spin states of 6 Li, \F,i71f) — 1 1/2, 1/2} , 
1 1/2, — 1/2), calculated from updated potentials consis- 
tent with the recently measured zero crossing of a at 
B ~ 530G |H El 

Beside the broad resonance near 
855 G on which we concentrate in this work, there is a 
very sharp resonance at 545 G, and increased atom losses 
have been measured peaking around 545 G and 680 G 
||. Recently, strong interactions in this Fermi mixture 
have been demonstrated at a fixed magnetic field of 930 G 
through anisotropic expansion of the gas [jj . In this Let- 
ter we study the atomic interactions through the entire 
range of the Feshbach resonance, from 600 G up to 1.3 kG. 
We measure three physical quantities : the interaction en- 
ergy with positive and negative values, the anisotropy of 
the atomic cloud during the expansion, and the atom loss. 
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FIG. 1: Predicted scattering length vs magnetic field in Li 
\F,m F ) = 1 1/2, 1/2), |l/2, —1/2) mixture. 

We find the resonance at 800(40) G, where the trapped 
gas is most strongly interacting, with a ratio between the 
interaction energy and the kinetic energy reaching —0.3. 
Surprisingly, we observe a large shift (~ 80 G) between 
this position, and both the location of maximum loss in 
the gas and the change of sign of the interaction energy. 
Our results are in partial agreement with the physical 
picture of a Feshbach resonance in the region of strong 
interactions but also reveal unexpected effects. Finally 
they point towards the best experimental conditions to 
search for the superfiuid transition in 6 Li. 

Our experimental approach to measure the interaction 
energy of a two component Fermi gas is based on the 
analysis of time of flight (TOF) images of atoms released 
from an anisotropic trap. The energy of the trapped gas, 
prepared at any value of the B field, is the sum of poten- 
tial, kinetic and interaction energies, E pot + Ekin + E^ nt . 
Switching off abruptly the trapping potential (E pot — > 0) , 
the gas is released and expands for a variable time be- 
fore an absorption image is recorded. At long expansion 
time, the spatial distribution of the cloud reflects the 
velocity distribution. This procedure is done routinely 
for BEC studies Q] and has been recently investigated 
theoretically for Fermi gases in ^3- The novelty of our 
approach is in what follows. Because of the low induc- 
tance of the coils used to produce the magnetic field B, 
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this field can be switched off rapidly (< 20 jus) at any 
desired time during the expansion of the atomic cloud. 
For B = the atoms have negligible interactions since 
a ~ (Fig.[TJ. As a consequence, the expansion of the 
gas can be recorded without atomic interactions during 
the time of flight period (B — 0), or with interactions 
during the TOF period (B ^ 0). Since expansions with 
B = are ballistic, they reflect the kinetic energy of 
the initially trapped gas, E^in- On the other hand, for 
TOF with B ^ 0, the interaction energy is converted 
into kinetic energy during the early stage of the expan- 
sion (< 150 /is). TOF images at long time reflect the 
released energy, E re \ = Eki n + Ei nt jig . Comparisons 
between TOF images with B = and B ^ allow us 
to simply deduce the ratio E int /Ekin- We operate at 
temperatures between 0.5 Tp and Tp, where Fermi de- 
generacy does not play an important role. Even in this 
nearly classical regime, the gas is found to be strongly 
interacting. The observed expansions with B ^ are 
anisotropic because of collisonal hydrodynamics, as pre- 
dicted in [HEI, and observed in [fl Ell Ejf . 

Our experimental setup has been described previously 
0, Hj|. A gas of 3 x 10 5 6 Li atoms is prepared in the 
absolute ground state |l/2, 1/2) in a Nd-YAG crossed 
beam optical dipole trap at a temperature of 10 /iK. The 
horizontal beam (resp. vertical) propagates along x (y) , 
has a maximum power of 2.5 W (2.5 W) and a waist of 
~ 25/itm (~ 30 fim). At full power the trap oscillation 
frequencies are u> x /2ir = 2.2(2) kHz, Lo y /2n = 3.0(3) kHz, 
and uj z /2tt = 3.7(4) kHz, as measured by parametric ex- 
citation. Using a radio frequency field, we drive the Zee- 
man transition between 1 1/2, 1/2) and 1 1/2 , —1/2) to pre- 
pare a balanced mixture of the two states at any chosen 
value of B between 2 G and 1.3 kG. Using a Stern-Gerlach 
method we check that the populations in both states are 
equal to within 10%. 

We observe that atom losses occur with very different 
rates below and above resonance. Near 720 G the life- 
time of the gas is on the order of 10 ms, whereas near 
900 G, the lifetime is in excess of 10 s. The latter is sur- 
prisingly large in comparison with similar situations for 
bosons near a Feshbach resonance [24| • Therefore we per- 
formed two sets of experiments : one with the spin mix- 
ture prepared above resonance (1060 G), and the other 
below resonance (5G). In the first set, evaporative cool- 
ing is performed by lowering the power of the vertical 
beam. It produces 2N = 7 x 10 4 atoms at T ~ 0.6 Tp 
with k^Tp — hui (6./V) 1 / 3 and u) — (ujxUjyUJz) 1 ^ 3 . The trap 
is nearly cigar shaped with frequencies uj x /2t: = 1.1 kHz, 
uiy/2ir = 3.0 kHz, and uj z /2tt — 3.2 kHz. The magnetic 
field is adiabatically ramped in 50 ms to a final value 
where TOF expansion images are taken. The effect of 
strong interactions is illustrated in Fig.EJ where both 
typical images and a plot of the expanded cloud gaus- 
sian sizes r x and r y are displayed. In Fig.^, expansions 
with B = reveal the initial momentum distribution of 



the cloud which is isotropic as expected. We deduce the 
total kinetic energy E^m of the trapped gas mixture from 
gaussian fits to this distribution. For our moderate quan- 
tum degeneracy, the temperature can be estimated from 
k B T = 2E km /3. We find T ~ 3.5 0.6 Tp, constant 
for fields between 0.8 kG and 1.3 kG. 
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FIG. 2: Left: Geometry of the trapped atomic cloud (center) 
and expansion images without magnetic field (a) (B = 0), 
and with magnetic field (b) (B 7^ 0). Right: Corresponding 
gaussian sizes of the expanded clouds along x (squares) and 
y (open circles) vs magnetic field, a : Time of flight images 
after 650 [is expansion with B = 0. b: Images after 400 fis 
expansion with B/0 and 250 [is with B — 0. Images shown 
on the left are for B = 900 G. 

In Fig-EI 3 ; expansions with B ^ are anisotropic. Lit- 
tle expansion is seen along the weak axis of the trap. The 
ellipticity of the cloud is inverted because of hydrody- 
namic behavior during the expansion • Collisions then 
redistribute the gas energy in the direction of maximum 
density gradient. The anisotropy r y /r x ranges from 1.1 
at large B field to 1.4 near 0.8 kG whereas the hydrody- 
namic scaling equations 19| predict an anisotropy of 1.53 
in the fully hydrodynamic regime. The usual criterion for 
hydrodynamicity is found from the ratio R of the mean 
free path Ao = {n^a)^ 1 , where uq is the peak density, 
over the radial size r ra d of the cloud. For a classical gas 
neglecting mean field interactions we find : 



R 



A (2tt) 3 / 2 (k B T\ uj rad 



■V) 
T*rad 



N<7 



R <C 1 (R 3> 1) corresponds to the hydrodynamic (col- 
lisionless) regime. With the predicted value of a = 
-185 nm at 1060 kG, and a = Ana 2 , R = 0.03. Therefore 
in the early stages of the expansion, the gas is hydrody- 
namic as in (jf |2l|, |22j ■ Furthermore, with these large val- 
ues of a and our typical temperatures, the scattering cross 
section is energy dependent. At B = 1060 G, \ka\ — 0.95, 

where k = \J rnksT '/2h 2 is the typical relative momen- 
tum of two colliding atoms. The cross section is reduced 
and becomes unitarity limited, a = Aira 2 /(l + k 2 a 2 ) ~ 
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in/k 2 for \ka\ 3> 1 26]. Consequently, as the magnetic 
field is decreased below 1.3 kG, the gas gradually enters 
deeper in the unitarity regime. 

This effect has an important consequence on the gas 
behavior during hydrodynamic expansion: since the rel- 
ative momentum k of the colliding atoms decreases as 
k oc n 1 / 3 (where n is the density) |25j, a increases. For a 
spherical cloud, the hydrodynamic factor R would remain 
constant during expansion while for a cigar shaped cloud, 
the expansion is mostly 2D, and R decreases like n 1 / 6 . 
The gas becomes more hydrodynamic as the expansion 
proceeds until this model breaks down when relative mo- 
menta k become too small to remain in the unitarity limit 
|27j . Then R increases as n -2 / 3 for a spherical geome- 
try, and as n^ 1 / 2 for a cigar. Consequently, the larger 
\a\, the longer the expansion remains hydrodynamic and 
the stronger the anisotropy is. When the B field is de- 
creased in Fig.0 the anisotropy increase indicates that 
the scattering length becomes more and more negative. 
Therefore, this puts an upper bound of ~ 800 G for the 
position of the Feshbach s-wave resonance. 

In the second set of experiments, we focused on the 
lower magnetic field values between 550 G and 820 G, 
the region where losses have been observed. The mixture 
is prepared at low magnetic field, and the field is ramped 
up to 320 G, where the scattering length vs B has a 
local minimum of a = —8 nm. Evaporative cooling is 
performed there, leading to T = 2.4 /iK and to uj x /2tt — 
0.78 kHz, ujy/i: = 2.1kHz, and lo z /2tt = 2.25 kHz. Fi- 
nally, the B field is ramped to 555 G in 50 ms where 
a ~ 0, and then in 10 ms to different values between 
600 G and 850 G where TOF expansions are recorded. 
The gaussian widths and number of detected atoms vs B 
are plotted in Fig.|3 Comparing with the data of Fig.0 
we observe several new features. First, the number of 
detected atoms has a pronounced minimum near 725 G, 
a position compatible with previously published results 
dealing with losses |5j]. These losses are associated with 
a strong heating of the cloud (Fig.EH). 

Second, in figure 0(b) for B < 675 G, the expansion 
is isotropic, consistent with a collisionless regime. Here 
a is predicted to be < a < 55 nm. Since ka < 0.35, 
the scattering cross section is independent of energy and 
the gas is not hydrodynamic (R > 1). Above 700 G, a 
pronounced asymmetry in TOF images appears abruptly. 
The gas enters the hydrodynamic regime (and unitarity 
limit) with the same inversion of ellipticity as in Fig.[21 a 
signature of strong interactions. We expect that on res- 
onance, the elastic cross section is unitarity limited for 
all relevant values of energy. Therefore, the maximum 
anisotropy, which we measure at B — 800(40) G, locates 
the peak of the Feshbach resonance, in agreement with 
the predicted position of 855(30) Gauss. A unique fea- 
ture of the resonance is the very large shift (~ 80 G) 
between the resonance peak and the maximum of loss 
and heating. This can be qualitatively explained by the 
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FIG. 3: a and b: Gaussian sizes along x and y vs magnetic 
field. The detection conditions and symbols are the same as 
in Fig.|5| a: Expansion with B = 0. b: Expansion with 
B c: Number of detected atoms vs magnetic field. Open 
triangles and black diamonds correspond to graph a and b, 
respectively. Line LI corresponds to the maximum of loss, 
and L2 to the resonance position. 

creation for a > of weakly bound molecules by three- 
body recombination. According to [28j], in this process 
the binding energy of the molecule Ti 2 /ma 2 for large a, 
is dissipated into kinetic energy of the atom + molecule 
system. For very large values of a, the binding energy is 
small and the heating associated with three-body recom- 
bination is negligible, and the molecules remain trapped. 
For small or intermediate values of a, the collision prod- 
ucts are highly energetic, resulting in strong heating and 
loss. At the maximum of loss (720 G), a = 102 nm and 
h 2 /kBma 2 = 7 /iK is on the order of the trap depth in 
the weakest direction x. 

As explained above, the interaction energy can be cal- 
culated directly from the difference of release energies, 
obtained from the gaussian sizes of Figs.[21and|3| In Fig.0] 
is plotted the ratio of interaction energy and kinetic en- 
ergy vs magnetic field. As B is increased above 550 G, 
Ei n t first goes up, due to increasing a, in accordance 
with mean-field theory. Near 720 G, Ei nt changes sign 
abruptly. The ratio Ei nt /Eki n then exhibits a plateau 
near —0.12 up to 850 G. Note that, for the data at low 
field (Fig.|2J), the number of atoms and the temperature 
are not constant as a function of B. From the data at 
higher fields, we observe in Fig.0] that the interaction 
energy becomes less and less negative above 800 G. The 
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tween molecular condensates and superfluid Fermi gases 
are promising. 
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FIG. 4: Ratio of interaction energy over kinetic energy 
vs magnetic field. Open squares: atomic samples prepared 
above resonance as in Fig.[3] Crosses: average of three sets of 
data, recorded in conditions of Fig.|3] The difference between 
crosses and open squares in the overlap region is due to differ- 
ent experimental conditions described in the text. Solid line : 
Mean field theory calculated with a total number of atoms of 
7 x 10 4 and a kinetic energy of Ekin/ks = 5.25 fiK, values 
corresponding to the open squares conditions. The gray area 
indicates the region of losses. 

apparent discrepancy between the two sets of data in 
the overlap region can be explained by different experi- 
mental conditions and, in particular, different atom num- 
bers, temperature, and confinement. The line in Fig.^l 
shows the result of a mean field calculation without any 
fit parameters based on a full energy dependent scatter- 
ing phase shift. The mean field energy Emf is then ob- 
tained from a self-consistent distribution function that 
contains the self-energy and trapping potential. The 
kinetic energy is calculated from the same distribution 
function. Theory and experiment agree quantitatively for 
B > 850 G where a < 0, and for B < 720G, where a > 0. 
However, the observed shift between the predicted posi- 
tion of the resonance and the change of sign of the mean 
field 720(20) G is not reproduced by this theory. A pos- 
sible explanation could be that dimer molecules, which 
are expected to be present for a > and close to reso- 
nance, are not part of the calculation. A large negative 
partial mean field contribution due to atom-molecule and 
molecule-molecule interactions could shift the position of 
the transition from positive to negative mean field. 

In summary, we have studied a Fermi gas mixture in 
the strongly interacting regime near a Feshbach reso- 
nance. The position of the resonance is found in agree- 
ment with theoretical expectations. New features have 
been observed which may be the signature of richer 
physics, for instance molecule formation |29l| . Anisotropic 
expansions are observed both for repulsive and attractive 
mean field interactions, in a moderately degenerate Fermi 
gas. This is interpreted in terms of collisional hydrody- 
namics without invoking Fermi superfluidity. Prospects 
for producing deeply degenerate Fermi mixtures in the 
superfluid state and for investigating the transition be- 
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